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A mutual-Chern-Simons Lagrangian is derived as a minimal field theory description of the phase- 
string model for doped antiferromagnets. Such an effective Lagrangian is shown to retain the full 
1^-^ ■ symmetries of parity, time-reversal, and global SU(2) spin rotation, in contrast to conventional 

f*"*) ' Chern-Simons theories where first two symmetries are usually broken. Two ordered phases, i.e., 

, antiferromagnetic and superconducting states, are found at low temperatures as characterized by 

. "dual" Meissner effects and dual flux quantization conditions due to the mutual-Chern-Simons gauge 

structure. A "dual" confinement in charge/spin degrees of freedom occurs such that no true spin- 
charge separation is present in these ordered phases, but the spin-charge separation/deconfinement 
serves as a driving force in the unconventional phase transitions of these ordered states to disordered 
states. 
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Gauge theory description has become essential in studying doped Mott insulators. The physical necessity may be 
traced to the Hilbert space restriction in a doped Mott insulator. For instance, the high-T c cuprate superconductors 
at half-filling are believed to be an antiferromagnetic (AF) Mott insulator 1 , in which the charge sector at low energy 
is totally frozen up by the Coulomb interaction. After doping, the low-energy charge degrees of freedom do emerge, 
but remain highly restricted in the Hilbert spacei. To characterize such a Hilbert space restriction, a spin-charge 
separation description, namely, by introducing 2,3,4 spinless "holon" of charge +e and neutral spin-1/2 "spinon" as 
the essential building blocks of the restricted Hilbert space, has become an effective and useful way. Here "holons" 
£^ ' and "spinons" do not necessarily turn out to be true low-lying elementary excitations in the end, because generally 
local gauge field(s) will emerge^ to mediate interactions between these "holons" and "spinons" , and may even lead 
7-H ' to the confinement of them if either a true spin-charge separation does not exist or the decomposition is not done in 
Q\ ' a correct way. In general, one always ends up with a gauge theory description for doped Mott insulators where the 
gauge interaction can greatly influence the low-energy dynamics of the charge and spin degrees of freedom. 

Several kinds of (2+l)-dimensional gauge theories have been proposed for doped two-dimensional (2D) spin-1/2 
antiferromagnets related to the high-T c cuprates. A U(l) gauge theorjiSii based on the slave-boson approach to the 
t — J model is one of the most intensively studied. Its gauge structure may be directly visualized by noting the gauge 
invariance of the electron operator in the slave-boson decomposition^ 



bin, (i) 



. under a U(l) transformation: bi — > e l8i bi and — ► e iei fi a , where bi denotes the bosonic "holon" operator and fi a 
the fermionic "spinon" operator. Along the same line, the SU(2) non-Abelian gauge theories^ and Z2 gauge theories^ 
O ■ have also been proposed and studied. 

O ■ The slave-boson approach is considered to be convenient in dealing with the superconducting (SC) regime but has 
less advantage in describing the AF state near half-filling. On the other hand, gauge theoriesi2*ii*i2ii2ii2ii£ based 
on the slave-fermion, Schwinger boson decomposition are believed to be useful in studying a lightly doped AF state. 
Here the electron operator is written as^ 



C la = ftb ia (2) 

where /$ denotes the fermionic "holon" operator and bi a the bosonic "spinon" operator. Besides the slave-boson 
and slave-fermion decompositions, slave-anyon decompositions have also been investigatediiiSii^iSS. Different gauge 
structures mentioned above originate from different decompositions and/or different mean-field decouplings. But a 
common feature for these gauge theories is that both "holon" and "spinon" share the same gauge field. 

Recently, a different gauge theory description has been constructed 21 based on a distinctive decomposition of the 
electron operatos 2 ^^ 
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which is known as the bosonizationS^ or phase string decomposition^ because holon and spinon operators, h\ and 
6j CT , are both bosonic, with the fermionic commutations relations of the electron operator being restored by the phase 

string operator, e l ® icr = (— cr^e'st** - "^* L Here internal gauge invariance appears as U(l)xU(l): hi — > e^hi and 
$J — > $^ + 2<^; 6 l(T — > e 4CTX »& i(T and $^ — > $^ +x<- Consequently there exist a pair of U(l)xU(l) gauge fields coupling 
to the holon and spinon fields, respectively, in the resulting gauge theory, called the phase string model, derived^! 
based on the decomposition © and the bosonic resonating-valence-bond (RVB) mean-field saddle-point, where the 
normal U(l) gauge freedom 24 (like the one in the slave-boson case) is broken by the mean-field decoupling 25 . 

In the slave-boson (or slave-fermion) U(l) gauge theory, the external U(l) gauge field (i.e., the electromagnetic 
field) couples to both holons and spinons&i, thanks to the same internal U(l) gauge field they share. So both holon 
and spinon carry some fractions of the electron charg o 7 ! 26 . In contrast, in the phase string model, the external 
electromagnetic field only couples to the holon degrees of freedom, without being directly transferred to the spinon 
part as the latter sees a different gauge field. In this sense, the holon carries the full charge of +e in the phase string 
model. 

Without a bare kinetic energy, the single U(l) gauge field in the slave-boson (or slave-fermion) theory fluctuations 
strongljii^, which makes the theory a strong-coupling one. On the other hand, the U(l)xU(l) gauge fields are 
topological ones with their strengths constrained to the densities of two matter fields (see Sec. II) such that their 
fluctuations are much more mildly, suitable for a perturbative treatment. In particular, the no-double-occupancy 
constraint of the doped Mott insulator, which is enforced by the violent gauge fluctuations in the slave-boson (or 
slave-fermion) theory, is realized in the phase string model in a quite different way. Namely, the U(l)xU(l) topological 
gauge fields will introduce mutual repulsions between holons and spinons, where holons perceive spinons as vortices 
and vice versa. As it is well known, a particle cannot go to the core of a vortex of its own field where the density 
of such a matter field vanishes. In the phase string model such a vortex core of one species is always occupied by a 
different species such that the no-double-occupancy is naturally enforced. 

Furthermore, the weak (logarithmic) confinement of spinons and holons at low energies and low temperatures has 
been also founoj 2 ^ 2 ^ 2 ^ in the phase string model, as opposed to the strong confinement in usual 2D compact U(l) 
gauge models in slave-boson or slave-fermion theory^S^i. In the latter, an effective gauge theory may have a serious 
infrared divergence^ which makes the gauge theory very difficult to deal with mathematically. The former is usually 
much more manageable than the latter in this regard. 

However, the Hamiltonian formalisms! of the phase string model, in which a gauge field seen by one species is 
constrained to the density (number) of different species, is not very convenient for studies beyond the mean-field level. 
In this paper, we shall develop a Lagrangian (path-integral) formalism of the phase string model. We show that the 
effective low-energy Lagrangian describes two matter fields, holon and spinon, minimally couple to two different U(l) 
gauge fields. These gauge fields have no their own kinetic terms either, but there is a mutual- Chern- Simons term 
which entangles two gauge fields together. We call this as a mutual-Chern-Simons description, which constitutes a 
minimal field-theory description for the phase string theory. 

The gauge structure of such a (2+l)-dimcnsional mutual-Chern-Simons theory is very unique in many aspects as 
compared to the gauge theories proposed before. We demonstrate that the physical symmetries, which include parity, 
time-reversal, and spin rotational symmetries, are precisely preserved in such an effective theory. By contrast, in 
usual Chern-Simons (anyon) theories^2i2i^4, the parity and time-reversal symmetries are explicitly broken, including 
the mutual-Chern-Simons theory previously proposed-^ for describing the double-layer quantum Hall effect system. 

We further show that there exist two low-temperature phases in such a theory at low doping. One is an AF state 
which recovers the AF long range order (AFLRO) of the Heisenberg model at half-filling and may survive at small 
doping concentration. The other is an SC state. Two phases are characterized by dual Meissner effects and dual 
flux quantization conditions, accompanied by a dual confinement, which are the direct consequences of the mutual- 
Chern-Simons gauge fields interacting with two matter fields when one of them experiences Bose condensation. Such 
a mutual duality connecting the AF and SC states or spin and charge degrees of freedom, is quite different from the 
usual duality descriptions proposed^S^LSi^ for the cuprate superconductors, where the conventional boson-vortex 
duality is used to describe an ordered-disordered transition. 

In the SC phase, for example, the Meissner effect and hc/2e flux quantization are similar to the predictions by a 
conventional superconductivity theory, and the spinons are found to be confined such that to drop out of the physical 
spectrum. Only integer spin excitations, as composed of confined spinon pairs, are allowed in the bulk state. But 
as a unique prediction, a single spinon (an S = 1/2 moment) does appear in the center of a magnetic vortex core. 
It forecasts that the spin fractionalization will occur in the pseudogap phase, as the latter may be viewed as the 
proliferation of the vortex core state above the superconducting transition T^Z^. 

In the AF phase, on the other hand, the spinon condensation may be viewed as a two-component "superfluidity" . 
The dual Meissner effect means that a holon is an "alien" object in the spinon condensate, and the dual flux quanti- 
zation condition corresponds to that a meron (vortex) is produced in the spinon condensate to which a holon must be 
confined to, just like a spinon is confined to a magnetic vortex core in the above-mentioned SC state. As a result, only 
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the "neutral" object of a holon-meron composite, not the holon itself, appears in the low-energy physical spectrum, 
which has a dipolar spin configuration at long distance, coexisting with the AFLRO in a dilute hole concentration 
regime. 

The remainder of the paper is organized as follows. In Sec. II, we briefly introduce the effective Hamiltonian of the 
phase string model. In Sec. Ill, we first derive the Lagrangian (path-integral) formalism in the lattice version. Then we 
obtain the low-energy mutual-Chcrn-Simons gauge theory description in the continuum limit. In Sec. IV, we examine 
the symmetries, including parity, time-reversal, and spin rotational symmetries, of the mutual-Chern-Simons theory. 
In Sec. V, we study two low-temperature ordered phases based on the mutual-Chern-Simons theory and discuss 
how holons and spinons behave in the AF and SC phases, respectively, where dual confinement of holons/spinons is 
revealed. Finally, the conclusions are given in Sec. VI. 



II. PHASE STRING THEORY: A MINIMAL MODEL OF DOPED ANTIFERROMAGNETS 

The phase string theory has been proposed 21 '— as a low-energy effective description of the doped antiferromagnets 
at low doping. The "minimal" Hamiltonian of the phase string theory is composed of two terms, i? s tring = Hh + H s , 
in which the charge degrees of freedom are characterized by the "holon" term 



H h = -t h J2 ( elAh ) h \ h o + H.c. (4) 

m 

where th ~ t and the "holon" operator, h\, is bosonic; The spin degrees of freedom as described by the "spinon" term 

H s = -J s J2{e laA -)blbU+H.c. (5) 

where J s ~ J and the "spinon" operator, b\ a , is also bosonic. Here the gauge fields AL and are decided by the 
topological constraints: 



C i6S c 
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where n h liy and denote the "spinon" (with index a) and "holon" number operators at site I, respectively, and the 
path C is an arbitrary loop made of the nearest-neighbor (nn) links with Sc denoting the area enclosed by C. 

Basic features of this model are as follows. At half filling, the gauge field can be set to zero in (JSJ and H s 
reduces to the Schwinger-boson mean-field Hamiltoniani^, which describes both the long-range and short-range AF 
correlations fairly well. Upon doping, A^ is no longer trivial due to constraint ©, which describes that each "holon" 
behaves like a 7r-fluxoid as felt by the "spinons". Thus, A^ will play the role of dynamic frustrations, introduced 
by doped holes, that acts on the spin degrees of freedom. Similarly, the "holons" are also subjected to dynamic 
frustrations, from the spin background, via the gauge field Afj in @. The spin and charge degrees of freedom are 
thus mutually frustrated in the phase-string model in terms of two topological gauge fields, A^ and Ajj . 

The phase string model outlined above incorporates, as a minimal model, three most essential characteristics of 
the doped antiferromagnets described by the t — J model. They are: (i) the restricted Hilbert space of doped 
Mott insulators, which is characterized by the spin-charge separation formalism with holons and spinons as basic 
building blocks; (ii) strong short-range AF correlations as provided by the bosonic RVB description in (JSJ, which can 
naturally grow into an AFLRO state as the doping concentration is reduced to zero; (iii) the mutual singular influence 
between the charge and spin degrees of freedom as represented by two topological gauge fields, A^ and AL, which 
mathematically capture the phase string effect identified 2 ^ in the t — J model. Such a mutual interaction has been 
s } 10wn 2i i 27 i 28 4 29 ^ Q responsible for some nontrivial physical properties of the model in close connection with the 
high-T c materials. 

In the phase string formalism, the spin operators are expressed in terms of the spinon operators in the following 
nontrivial forrn^ 
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FIG. 1: A regularization of the contraints in © by introducing dual lattices is shown. A spinon (denoted by an arrow) and a 
holon (denoted by an open circle) stay in dual lattices (solid and dashed ones, respectively) , with the gauge fields Ah and Aj j 
defined on the links of two dual lattices, respectively. The closed loop, C (C*), of a spinon (holon) can be arbitary without 
crossing holons (spinons). See text for the detail. 



St = (Srf = (-l)%b tl e^ (7) 

where the phase <£>f appearing in Sf can be decided by the relation <£>' 1 — <£>^ = 2 Ay- for two nn sites, i and j, which are 
not occupied by the holes. Under this definition, the spin operators as well as the effective Hamiltonian are invariant 
under the gauge transformation b ia — > bjo-e" 7 ^* , <F' 1 — > <J>^ + 2(f> i: Ay — > Ay- + <f) i — <j>^. The holon-dependent phase 
factor $^ in Q further illustrates the intrinsic mutual entanglement between spin and charge degrees of freedom in 
the phase string theory. 

III. MUTUAL-CHERN-SIMONS GAUGE-THEORY DESCRIPTION 
A. Lagrangian Formulation 

The treatment of the Hamiltonian formalism of the phase string model may not be convenient beyond the mean- 
field approximation because the gauge fields, Ay and Ay, defined in ©, are themselves operators depending on 
the dynamics of the matter fields. In order to deal with the phase string model [Q and JSJ] more conveniently, a 
Lagrangian (path-integral) formalism will be introduced in this section. 

First of all, let us re-express the topological constraint @ locally. As pointed out before, the original no-double- 
occupancy constraint in the t — J model can be realized in the phase string model by the mutual repulsion between 
spinons and holons via A^ and Ay. As a consequence, the closed path C of a holon/spinon in © will not cross 
spinons/holons and thus effectively avoid a singularity occurring when a spinon and a holon simultaneously stay at 
the same site (as each spinon/holon carries a 7r-fluxoid seen by a holon/spinon). Following this, then, it is physically 
reasonable to implement a regularization in the topological constraint (jSJ) by introducing two sets of dual square 
lattices, respectively, for spinons and holons to stay, as illustrated in Fig. ^ In this way, a closed path C for spinon 
and a C* for holon on different lattices can be arbitrary without worrying to cross the opposite species, either holons 
or spinons. Presumably no important low-energy physics will get lost by such a local regularization. 

Here and below, the minuscule (majuscule) Latin letters i,j (I,J) will be used to label the dual lattice sites for 
spinons (holons). The Greek letters a, (3, 7 will be used for 2D spatial indices, 1 and 2, while /i, v, A for the three 
dimensional space-time indices, 0, 1, 2. Then the topological constraint © can be re-expressed in a compact form as 
follows: 

e af3 A a A h p (i) = nn 1 } 
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(8) 



in which the link fields A^(i) = A^ +Sl i and A s a (I) = Aj T _ & , with a = x, y, and the difference operators, A a , on the 
two sets of the dual lattices are defined by A a f(i) — fit + a) — f(i) and A a f(T) = f(I) — f{I — a), respectively. Note 
the slightly different definitions of link variables and lattice difference operators on two dual lattices, so as to keep 
the symmetric forms in (JSj. 

In the path-integral formulation, the topological constraint (jSJ can be enforced by introducing two Lagrangian 
multipliers, A^(i) and Aq(J), as follows 



n h _ - e <*PA a A h p(i) 



E 



cm. 



-e a/3 A a A%(I) 



(9) 



Once the topological constraint is implemented by the Lagrangian multipliers, the gauge fields, A s Ia and A^ a , can be 
treated as independent gauge variables in the Lagrangian formalism. In order to get the correct form of the Lagrangian 
for this system, we need to first identify the canonical momenta of the gauge fields, A'^ and Aj j. 
It is helpful to consider the continuity equation for the holon density: 



d t n'} + A a 4 +&I = 0. 



Using the topological constraint in 



and the definition of the conserved holon current, Jj +& j 



(10) 



-T-p — one gets 



-e^A a A h Ji) 



6H< 



string 



SAl 



such that [under a proper gauge choice of ^4«(i)] 



dtA h (i) 



string 



S^TT-^A^I))' 



(11) 



(12) 



Equation (|12|l is just the canonical equation of motion for A^(i), and one can thus identify the canonical momentum 
Il£(i) = --e^A^ (I). In other words, the spatial components of the gauge fields A h and A s are canonically conjugate 
to each other. (The temporal components, Aq and Aq , have no canonical momenta since they do not have independent 
dynamics in the above formulation). 

Following the standard canonical quantization procedure, the Euclidean Lagrangian (with the Wick rotation t — > 
—it) of this system can be derived straightforwardly as follows 



Lstn„ g = En^HH^M + E^c. 6 - 

i i 

= L h + L s + Lcs 



E h \ d ohi 



H. 



string 



where 



Lh = 



L s = 



E h\ [do - iA s (I)] hj - t h E (hW A '-'hj + h. 
i (U) 



E bl [do laAlii)} b ia - J s E (ble^b]_ a + h.c. ) 



^E h \ h i - NS 
a(E 6 1 6 « 

\ icr 



N(l-8) 



Lcs = ^e^A^I^A^) 



(13) 



with do = d T . Note that one can also use a procedure similar to p0(l - l|12(l to define a conjugate field 11^(7) — 
— ^e a PAp(i) for A s a (I) and the resulting Lagrangian remains the same as above. 

Therefore, the Lagrangian formalism of the phase string model describes that the two matter fields, bosonic spinons 
and holons, are minimally coupled to U(l) x U(l) gauge fields, A 8 ^ and A 1 ^, whose gauge structure is decided by the 
mutual-Chern-Simons term Lcs i n <|13|1 . In the following, we shall further derive the long- wavelength, low-energy 
effective Lagrangian based on such a lattice model. 
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B. Low-energy effective theory 

The Lagrangian (|13f) is written in a lattice form. It can be further simplified and reduced to a continuum version 
in the long- wavelength, low-energy limit. The procedure given below is quite standard and straightforward. 

Let us first consider the spinon Lagrangian L s , in which some careful treatment is needed in taking the continuum 
limit. We shall derive its low-energy action in the CP(\) formalism 41 by integrating out the short-range ferromagnetic 
fluctuations. 

First of all, we divide the square lattice into two sublattices, A and B, and redefine the spinon operator bi a at B 
sublattice as bi a . Then L s in (|13|l can be rewritten as 

L s = Y, bt(d -io-A%)b ia + Y bl(d -iaA%)b i<7 
J, H'Le^bL^ + h.' 

i^A,j— nn(i) ,a 




-N(l-5)\. (14) 

yi£A,a iEB,cr 

As usual, we introduce the following continuum fields 

b ia = z a (vi) + ir a (r t ) 
b\ +fl _ ff = z a {r l + fja) - 7r CT (ri + r)a) (15) 

in which i G A, fj = x, y, and a is the lattice constant. Then, by expressing the Lagrangian (|14J) in terms of z a and 
7T (j and taking the continuum limit a — > with A^(i) — > aA^(r), we obtain L s — J d 2 r£ s , in which 

C s = [js\(d a -iaA h a )z a \ 2 + a- 2 {\-U s )\z <y \ 2 ~\ [ - Aa~ 2 (1 - <5) 

a 

+ [- J « I ( 9 « - i(jAh o) ^-| 2 + fl_2 ( A + AJ s) 

a 

+a- 2 Y [K (do ~ «rA%) z a - n a (do + iaA' l ) z* a ] (16) 

(7 

By further integrating out the high-energy field 7r CT , we arrive at 



(17) 



Define the spin- wave velocity c s = ^/J S (A + 4J s )a and redefine the temporal components: r — > c s xq, Aq — > c s Aq, 
the low-energy effective action for the spinons can be finally written as 

S s = J d 2 rj^ ^ dx ^ -iaA*)z a \ 2 + m 2 \z a \ 2 ] , (18) 

Here the summations over /j, = 0, 1, 2 and a =|, J. are omitted for simplicity and the constant term Xa~ 2 (1 — 6) is also 



dropped. The coupling constant g = c s /2J s (1 — 5) , and the mass m s — c s 1 \J A 2 — 16 J 2 , in which A is decided by the 

spinon number constraint J d 2 xJ2a \ z &\ 2 = Na 2 . Note that here the z a field has been rescaled in the last step such 
that ^2 a \z a \ 2 keeps to be 1 per site on average even at finite doping. Therefore, in its final form, the long-wavelength 
theory for spinons consists of a massive, spin-1/2, and relativistic bosonic z a (spinon) coupled to a U(l) gauge field 

Ah 

The continuum versions of Lh and Lqs can be more straightforwardly obtained by directly taking the continuum 
limit a — > 0, with A 8 a (I) — » aA 8 a (r), Aq(I) — > c s Ag(r), and hi — > ah(r). The final form of the partition function can 
be written in the compact form 



Z = J DhDz 1 Dz l DA s DA h cxp I- J dx J d 2 r£ cS 
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in which 

£ cff = C h + C S + C CS (19) 

with 

C h = ftt [flb _ i{ Al + eAg)] ft + ftt ( ~^« ~ ^ ~ ft 

£ s = ^[l(^-i^KI 2 + Tn 2W 2 ] 

£ cs = V^A*^ (20) 

where in the holon Lagrangian density Ch, m-h — (2^a 2 ) -1 , is the vector potential of the external electromagnetic 
field, and — e is the electron electric charge. Note that the chemical potential /i in Ch has been absorbed into for 
simplicity. 

The Lagrangians in (|20|l constitute our final low-energy effective theory. They describe two matter fields, holons 
and spinons, minimally coupled to a pair of U(l)xU(l) gauge fields, and A 1 ^. The latter do not have their own 
kinetic energies, but are mutually "entangled" by the mutual-Chern-Simons term Ccs- Such a mutual-Chern-Simons 
term has been previously proposed— for describing the double-layer quantum Hall effect system. But here due to 
the fact that A h ^ couples to up/down spins with opposite "charges" in C S) the parity and time-reversal symmetries 
are explicitly retained (see below). The external electromagnetic field, A e , only directly couples to the holon field, 
indicating that the latter is the primary charge carrier (consistent with the definition of the holon) . This is in contrast 
to the usual U(l) gauge theory based on the slave-boson approach, in which both holon and spinon share the external 
electromagnetic field as if each of them carriers a fractional part of the charge e (as the result that both of them see 
the same internal U(l) gauge field). 

In the following section, we shall carefully examine the symmetries of this effective Lagrangian with a particular 
attention to the parity, time-reversal, and spin SU(2) rotational symmetries. 



IV. SYMMETRIES 



The symmetries of the present mutual-Chern-Simons Lagrangian will be studied in this section. The following 
discussions will be based on the low-energy effective Lagrangian (|20() , although all of them can be easily generalized 
to the lattice formalism in 

First of all, we note that the U(l) c har g e X U(l)s z gauge invariance of S e s is obvious according to l|2U|). Consequently, 
the global U(l) c h a rge invariance of the holons ensures the conservation of the electromagnetic charge in this system. 
Also straightforward is the translational invariance in (2+l)-dimcnsions. In the following, we shall mainly focus on 
the parity, time-reversal, and spin rotational symmetries, and show that they are explicitly retained in the present 
mutual-Chern-Simons gauge theory, in contrast to ordinary Chern-Simons theories in which the parity and time- 
reversal symmetries are usually broken. 



A. Parity 



In (2+l)-dimensions, the parity transformation is defined as a reflection with regard to a spatial axis, e.g., 

x->-x, y -> y, t > t (21) 

It is straightforward to verify that the effective Lagrangians, Ch, C s , and Ccs, remain invariant, respectively, under 
the parity transformation 1)21(1 . if the matter fields and gauge fields transform under (|21|) as follows 



-a, 



h^h 



Ah _ Ah Ah Ah Ah _ Ah 

A s — i A s A s — > —4 s A s — » 4 s (99) 

The parity transformations of the fields in (|22|l can be determined as follows. For example, according to the property 
of angular momenta, a spin should transform as an axial vector, namely, S x S x , S y — > — S v , S z — > —S z under the 
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parity transformation (|21|) . Thus the transformation of the CP{\) field z„ should be z a — » Z- a . On the other hand, 
the gauge field A h transforms as an axial vector and A s as a polar vector in (|22[1 . Indeed, in order to keep the 
invariance of Ch and C s , the parity of ^4* and A 1 ^ should be identical to the charge current ifj = — SCh/SA^ and spin 
current j" = —SCg/SA^, respectively. Furthermore, the parity transformations o f A s u and A 1 ^ are also consistent with 
the classical equations of motion for the Chern-Simons fields obtained based on (|20|l : 

f = L^ dv A' x , j* = -e^d v A\. (23) 

The parity invariance of the mutual Chern-Simons term is also related to the fact that the gauge field A h transforms 
as an axial vector and A s as a polar vector, in contrast to an ordinary U(l) Chern-Simons theory. 



B. Time-reversal 

Under the time-reversal transformation, 

t -> -t, r a ^r ai (24) 

the Z& and h field will transform as usual spinor and scalar fields, respectively. Using the same procedure as given 
above in the parity transformation, we can determine 

z a — > (JZ*_ a , h —>■ h* 

Ah _ Ah Ah Ah 

~ * ^Oi A s a — > —A^ (25) 

under the time-reversal transformation (|24p. It can be then easily checked that the Lagrangian C e g = Ch + C s + Ccs 
is also invariant under the time-reversal transformation. 

The parity and time-reversal invariances of the mutual Chern-Simons Lagrangian l|2(J|l are in sharp contrast to the 
violations of both, separately, in an ordinary U(l) Chern-Simons theory. As noted above, A h as an axial vector and 
A s as a polar vector in the mutual U(l) x U(l) Chern-Simons theory are the key for the restoration of the symmetries. 
Note that the charge conjugate symmetry is meaningless here since the holon Lagrangian Ch is non-relativistic and 
anti-holons are not well defined. 



C. Spin SU(2) rotation 

The demonstration of the global spin SU(2) symmetry in the present formulation is less straightforward than the 
other symmetries discussed above. The underlying reason is that the spin operators are expressed in an unconventional 
way in terms of the hi = (btf, bi[) T doublet according to JZJ. 

Let us consider a global SU(2) spin rotation defined by U — exp(i0 • S) . In terms of Q, one finds U^biU — 
{a- A ) l e- ia ^ , l 7V e '°'/V CT3 ** ' /2 {a 3 )%. Correspondingly, according to the definition of the CP {I) fields in (JSJl, the 
doublet z = (zj, z\) T under the SU(2) rotation U is given by 

C/-^(r, T)U = e ^3* h ( r ,r)/2 e ^.e/2 e - i a 3 *"(r,r)/2 z(r) ^ (2g) 

in which 

^"(r.rj^^fr.r). (27) 

(Note that in the Hamiltonian formalism, the single-valueness of $f in the spin operators Q is ensured by the 
topological constraint on A^ according to iJBJ . In the path-integral formalism, Q h is determined by (|27|l , and we show 
in Appendix A that to have a finite contribution to the partition function, $ h must still satisfy the single-valueness 
constraint: A^lc = § c d^ h dx^ = § c 2A l ^dx fl — 2nir, with neZ for an arbitrary loop C.) 
The spinon Lagrangian C s can be rewritten as 

C s = ^-{(D^ D^z + mtzh} (28) 
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in which D^z = (9 M — ia^A'^j z = d^^e wr3 *'*/ 2 ,zY Under the transformation l(26[). D^z transforms as 

U- 1 (D^z) U = e^-WDpZ. 

Namely, D^z transforms as the basic representation of the SU(2) group, and the SU(2) invariance of the Lagrangian 
(128(1 is proved. Independent of z(r, r), Lagrangians Ch and Ccs are obviously invariant. Therefore, the global spin 
SU(2) symmetry is indeed preserved in the present mutual-Chern-Simons theory. 

V. TWO ORDERED PHASES AT LOW TEMPERATURES 
A. AF phase at low doping 

1, Half-filling 

Let us first consider Lagrangian 1)19(1 at half filling. Without the presence of holons, one can find A 1 ^ = and £ e ff 
reduces to a CP(1) model 

C cS ^C s = ^-{\d^z\ 2 + m 2 s \z\ 2 }. (29) 
^9 

The saddle-point solution of (|2"9*)l can be obtained by a standard procedure after integrating out the CP(1) z-field 
and then minimizing the resulting action with regard to m 2 s (here the constant term —m 2 /2g previously dropped in 
C s has to be included) as follow a 42 i 43 i 44 



d 2 k 1 

T^T^ 5 , = 1 (30) 

47T" 2 k z + ojz + mi 



where cu n — 2irnT, n — integers. With a proper regularizationi 3 . in 1(31) (I . the mass gap m s can be determined at small 
T as 

m s w Texp(-^-). (31) 
T 9 

in the so-called renormalized classical region, where i = - — - > (here g c = ^ with A denoting a cutoff parameter 
in the regularization) . 

At T = 0, the mass gap m s = 0, and a Bose condensation takes place in the ground state with (z) =^0, corresponding 
to an AFLRO lying in the x-y plane: \Sf\ = (— 1)* < z\ >< zy >, which can be easily destroyed by thermal 
fluctuations at any finite temperatures as indicated by m s > according to 1(31(1 . 

The energy scale of the mass gap m s is always much smaller than the temperature, i.e., m s <C T, at T ^ i. 
Thus, io n = 2imT (n > 1) is usually much larger than the mass gap, which means that the quantum fluctuations will 
become negligible in a sufficiently long wavelength and low energy regime, where one may only consider the purely 
static (semiclassical) fluctuations. In the region of m s < k < c s (3, the effective Lagrangian of the CP(1) field will lose 
the Lorentz invariance and becomes 

C s *^\Vz\ 2 . (32) 

Such an effective Lagrangian can be also obtained in the renormalized classical region by using the 0(3) nonlinear a 
model 4 ^. 

2. Low doping 

In a sufficiently small concentration of holes, if the AFLRO or the Bose condensation of the CP(1) spinor fields 
persists, then the renormalized classical Lagrangian ((32|) remains applicable, which should be simply modified to 
couple to the gauge field A' 1 according to (|2"U|) as follows: 

£. = ^|(V-«r 3 A h )z| a . (33) 
^9 
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On the other hand, holons are coupled to in Ch, and two gauge fields are then entangled by the mutual-Chern- 
Simons term Ccs [ see <EUJ)] , which can be rewritten, up to a boundary term, as 



C cs = ~z Ah ■ ( ES x 2 ) + Z A o BS ( 34 ) 
field" strength. By integrating out A h , then, the spin dynamics will become entangled with the holon dynamics as 



i 

TT TT 

where we introduce E s = 9oA s — VAg as the "electric field" strength for A^ and B s 1 = V x A s • z as its "magnetic 
shown below. 

First of all, the integration over Aq will simply lead to B s = 0. In the following, one may then choose a proper 
gauge: A s — and E s = — S7Aq. Next, by using 

|(V - i<J 3 A h )z\ 2 = \Vz\ 2 + 2A h ■ v s + (A h f \~z\ 2 
with z = (z^, z*^j and v s = ^ (?Vz- Vz^z) , one has 



1 i . h ( 1 „ i 

-cs 



C s + C cs = t^IVzI 2 + A h • ( iv s - -E s xz 



under the constraint \z\ — 1, which, after integrating out A' 1 , arrives at 

=4 (' V5 ' 2 " |VS ' 2 ) + 2^ (ES)2 + i (E " X 2) ■ V "' (35) 
Finally, by introducing a unit vector n defined by 



and by using 



n = z^crz 



J|Vnr = |Vz| 2 -|v 



the low-energy effective Lagrangian reduces to 

C cff = i(Vn) 2 + ^(E- 1- + iXV, + £,, (*>) 

where 



= T-eoi/An • x <9 A n. 

47T 

This low-energy Lagrangian describes how the bosonic holons, via Ch, and spin twists, with topological charge 
density ICq, are coupled to a Maxwell gauge field A 8 with the "photon velocity" c — oo, that is, in the absence of |B S | . 
The only effect of such a non-relativistic gauge field is then to induce a 2D Coulomb interaction between two types 
of charged particles, including holons and spin twists characterized by ICq. Noting E s = — VAg and integrating out 
Aq in (|36|l . a potential term will emerge in the effective action as 

V = q 2 h J d 2 vd 2 v' In |r - r'| ( Ph + K s ) (r) ( Ph + £g) (r') (37) 

in which p h = fvh and q\ = tt 3 /g 2 , accompanied by a charge neutral constraint enforced in the thermodynamic limit 
on the low energy states of such a 2D Coulomb gas system, namely 



J d 2 rK(r) + /Cg(r)]=0. 



(38) 
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Thus, a holon has to be "confined" to a spin twist, satisfying 1 + J d 2 r^n ■ d x ii x d y h — 0, which leads to the 
quantization condition of the winding number of the unit vector {n(r)} in spin space as follows 

Q s = d 2 r^-n ■ d x n x d v ii 
J 4tt 

= 4 ^ 

Namely, each holon will be bound to a "meron" , which is a spin twist of the unit vector n whose winding number is 
half of that for a Skyrmion. 

According to the condition J3SJ|, one expects to find equal number of holons and (anti)merons at low temperatures, 
which are paired by the logarithmic- attractive interaction in l|37() . An unpaired holon or (anti) meron will cost a 
logarithmically divergent energy and thus is forbidden to appear. In other words, in the AF phase, a bare holon can 
not exist alone, but has to be always confined to a spin topological configuration (meron). Such an effect in the spin 
ordered phase is called the "holon confinement" . Note that a holon itself will also carry a spin vortex according to 
(|7|). the composite object formed by the holon- meron pair actually corresponds to a spin dipolar configuration in the 
real spin space, as previous identified in the phase string modeiS&Si. Since the (anti)meron is a semiclassical object 
without a coherent quantum dynamics, the dipole as a bound pair of a holon and a (anti)meron normally cannot 
move coherently either. That is, the holon will be self-trapped near the core of the meron in space and the translation 
symmetry is spontaneously broken. 

With the increase of doping, i.e., the number of holon-antimeron dipoles, one expects to see a screening effect on 
the confining potential V . It has been previously found that eventually a confmement-dcconfincment transition can 
take place beyond some critical doping concentration, where the screened 2D Coulomb interaction becomes short- 
ranged^^. Once the bosonic holons are free, they will experience a Bose condensation and the resulting phase is an 
SC state as to be discussed in the following section. In the SC phase, there exists a duality correspondence of the 
quantization condition (|39[) . which will ensures the flux quantization condition there. Correspondingly, (|39|l may be 
called a dual flux- quantization condition. 

Finally we remark that the hole self-trapping at low doping, discussed in the present work, is in contrast to a 
conventional picture for single hole moving in the AF background based on the numerical studies of the t — J model^. 
In the latter case, the doped hole is found to have finite spectral weight and a coherent dispersion with the bandwidth 
comparable to J. The discrepancy may arise from the small sample sizes in exact diagonalization calculations: The 
phase string effect, which leads to the mutual Chern-Simons gauge fields, starts to play the role of self-localization 
only when the sample sizes become larger than the localization length scales*. Further investigations from both 
analytic and numerical approaches to clarify this issue are needed. Possible experimental implications of self-trapping 
for lightly doped cuprate have been previously discussed in the phase string models. 



B. Meissner effect and spinon confinement in SC phase 

Now let us consider the other ordered phase with the Bose condensation of holons, (h) ^ 0, whose ground state is 
a superconducting one2£ with the Meissner effect and charge 2e minimal flux quantization as shown below. 
With (h) ^ 0, C h in (HU reduces to 

C h = i Ph (d <f> h — Aq) + -j^— (V^ - A s - A e ) 2 (40) 
2m h 

by writing h(r) — y /p^e l ^' h ^ r K The Chern-Simons term (|2U[) can be rewritten as 

£cs = --A s -(E' l xz) + -^B' l -z (41) 

by introducing the "electric" field = doA h — VAq and "magnetic" field B 71 = V x A h for the vector potential A h . 

Firstly, the "magnetic" field B h = B' 1 • z can be determined after integrating out Aq in the partition function and 
one obtains the condition 

B h = B h • z = n Ph (42) 

which is uniform and fixes the spatial component A h , such that E h = — VAq. Secondly, after integrating out A s , the 
resulting effective Lagrangian takes the following form 

C eS = £ s + {^-) |E h | 2 - iA h Q h (43) 
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in which 

Q h = -e 0vX d v (d x 4, h - A\) . 

Finally, we integrate out Aq in (|43|) . For our purpose, instead of using the continuous version (|20|l of C s , we shall 
use a simpler but more precise form of the term involving Aq based on the original L s defined in l|13[l . which reads 

C a = -iA%p.(T)+C. K = 0) 

in which p s (r) = p^(r) — P\(r) with pj(r)[p^(r)] denotes the density of up (down) spinons. Then, after integrating 
out Aq, one obtains the following effective action in (2+l)-dimensional Euclidean space 

S cff = J d 3 x^ [C a (A h =0)]+J dx Vsc 

where 

^sc = Q 2 S J d 2 rd 2 r' In |r - r'| (p s + Q h ) (r) (p. + Q h ) (r') (44) 

with q 2 = j^-- Similar to the case in the AF phase, in the thermodynamic limit, there is a charge neutral condition 
enforced as follows: 

" d 2 v[p s {v) + Q h {v)) 
A T - N x + (2N VOI - —) , (45) 



in which Aj (N±) is the total number of spin-up (spin-down) spins, A vor = w- J d 2 ve a P d a dpcf) h denotes the total 
number of 2n vortices in the holon field, and 4> e is the total external magnetic flux $ e = J d 2 Y<P vX d v A\. 
Since A-j-, Aj_, and A vor are all quantized to be integers, we find the minimal flux quantization condition 

la&J = y (46) 

where $o = 2tt (= hc/e in full units) is the flux quantum for a charge e system. 

Therefore, the external magnetic flux is not allowed to present in the bulk (i.e., the Meissner effect) unless it is 
quantized in multiples of half flux quanta given in Q46|l. In particular, for a magnetic flux quantized at the minimal 
half flux quantum $f n i n , there must be a spinon trapped near the vortex core according to l|45|) by noting that the 
"charge" 2iV vor of the vortices produced by holon field is always in units of 2ir (i.e., $o)- 

However, free spinons are not allowed in the bulk in the absence of the external magnetic flux. Indeed, for $ e = 0, 
the "charge neutral" condition (|45l) reduces to — 2Vj_ + 2N VOI = 0. As the result, a single spinon excitation, with 
S z = (JVf — iVjJ /2 = ±1/2, will violate the "charge neutral" condition, which in fact will cost a logarithmically 
divergent energy as each spinon behaves like a half vortex. Hence, in the superconducting state, the spinons-vortices 
must be always paired up (confined) in the bulk by the logarithmic force given in (|44|l . To be noted, not only the 
spinons with different spin indices (up/down) can pair up, those with the same spin indices (up/up and down/down) 
can also pair up to satisfy the charge neutral condition by involving a holon phase vortex with N vor ^ 0, e.g., 
Aror = Ij-^T = 2, N± — or A vor = —1,N^ = 0, Aj = 2. Since A vor does not appear in the rest of the action, 
these two excitations of S z = ±1 are energy degenerate with the state S z = 0, A vor = 0, to form S = 1 triplet 
spin excitations, consistent with the spin rotation symmetry generally demonstrated before. For the same reason, the 
single spinon bound to a magnetic vortex quantized at ^ should have a free moment with S z — ±1/2, because of 
the freedom introduced by A vor . 

Thus, the superconducting phase in the present mutual-Chern-Simons theory is characterized by the holon con- 
densation and spinon (logarithmic) confinement. We have seen that the fractionalization of spins (a single spinon) 
does not directly appear in the bulk low-lying excitation spectrum, but does show up in a magnetic vortex core^*. 
The deconfinement of spinon- vortex pairs will eventually occur at the superconducting transition temperature T^^. 
Finally, we point out that the symmetry of the superconducting order parameter, which is expressed in terms of the 
electron operator J3J), is d-wave like as discussed previously in Ref4£. 
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C. Mutual duality of two phases 

The doping effect and the interplay between charge and spin degrees of freedom are characterized by a mutual- 
Chern-Simons gauge structure in this model, as discussed in previous sections. The mutually dual characteristics of 
these two phases are summarized by the following table. 





AF 


SC 


Bose condensation 


< z >^ 


< h 


Coulomb gauge field 




Ah 


"charged" particle of 
Coulomb gauge field 


holon 


spinon 


external source of 
Coulomb gauge field 


meron 


magnetic flux 


"charge neutral" object 


holon-meron pair 


a. spinon pair 
b. magnetic flux + a spinon 


dual flux quantization 


lei = \ 


|d>e 1 = *a = M 

1 mini 2 2c. 


dual Meissner effect 


holon confinement 


a. spinon confinement 
b. spinon bound to magnetic flux 



We have shown that, at low doping, the spinon condensation leads to a spin AF order and forces a "confinement" 
on the holon part, making holons self-localized to ensure the AFLRO. On the other hand, at a higher doping, the 
condensation of bosonic holons forces a "confinement" on the spinon part, resulting an SC phase coherence. 

There are several distinctions between the two ordered phases. In the AF phase, the spinons condensate is a kind 
of two-component "superfluid" . Consequently the global symmetry is broken from SU(2) to U(l). In contrast, the 
ground state in the SC phase is a condensation of a scalar field - holons. As a result the global U(l) symmetry of the 
charge part is broken. 

Besides such a fundamental distinction, two phases share some common features originated from the duality in 
the mutual-Chern-Simons gauge structure. In both the AF and SC phases, there exist induced Maxwell terms that 
have only "electric field strengths" without the Lorentz invariance. There are "charges" coupling to these Coulomb 
gauge fields, including quantum particles (holons and spinons) and external sources without quantum dynamics. The 
"charge neutral" condition and the 2D Coulomb interaction among the "charged" objects lead to dual Meissner effects: 
In the SC case, an external magnetic flux as an external source must be quantized, and in order to realize a minimal 
quantum hc/2e, a single spinon must be bound to such a magnetic flux to form a "charge neutral" object; In the 
AFLRO state, as an external source, a (anti)meron is allowed as a topological excitation from the spinon condensate 
with a quantized winding number \Q S \ = \ and a holon cannot live alone and must be bound to such an external 
source to form a "charge neutral" object with a spin dipolar configuration. 

Finally we emphasize that the mutual-Chern-Simons theory in this work involves a mutual duality between the 
charge and spin degrees of freedom rather than a usual duality. A usual dual description has been also widely 
used2£i2Zi22ii22 in studying the doped Mott insulators, which deals with an ordered phase and the transition to a 
disordered phase in terms of the corresponding topological defects on the dual lattice. In a conventional dual-theory 
description, normally the AF and SC phases are not directly related. By contrast, in the mutual duality discussed in 
the present theory, the vortices of one species (holon/spinon) under condensation are themselves quantum objects of 
another species (spinon/holon), and two ordered phases, i.e., AF and SC states, can be naturally unified together—. 
A similar duality at low doping has been also investigated^ by starting from the slave-fermion approach— Based 
on the same phase string decomposition © but with a slightly different mean- field decoupling (see Ref M>) , a mutual 
duality between the AF and SC states has been also discussed recently in Refi^i. 

VI. CONCLUSION 

In this paper, we studied a new class of nontrivial (2+l)-dimensional gauge field structure - the mutual-Chern- 
Simons theory. The Lagrangian of such a mutual-Chern-Simons theory is derived as an effective low-energy description 
of the phase-string model for doped Mott insulators. This effective Lagrangian retains the full symmetries of parity, 
time-reversal, and global SU(2) spin rotation, in contrast to the conventional Chern-Simons theories where first two 
symmetries are usually broken. 
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The mutual-Chern-Simons theory as a minimal model for doped Mott insulators has a unique mutual duality 
structure. Two ordered phases found in this theory, the AF and SC states, are connected by dual Meissner/flux 
quantization effects and dual confinement/deconfinement. Namely, holons become vortices in the spinon condensed 
AF phase and spinons become vortices in the holon condensed SC state. The former leads to the holon confinement 
(a holon bound to a spin meron twist to form a "neutral" dipolar structure) and the latter leads to the spinon 
confinement and flux quantization (a spinon bound to a magnetic flux quantized at hc/2e to form a "neutral" object 
in the holon condensate). 

Such a mutual duality structure between the charge and spin degrees of freedom determines the essential competition 
between two degrees of freedom and provides driving forces for phase transitions to each other or to other disordered 
phases. Dual confinement means that there is no true spin-charge separation is present in these ordered phases since 
one species (spinon/holon) is always confined at low temperatures while the other (spinon/holon) is condensed. But 
the dual dcconfincmcnt will play an essential role in the transitions to disordered phases or at the boundary between 
two ordered phases where a quantum critical point may exist^S*^. The systematic evolution of the phase diagram at 
low doping is currently under investigation based on the mutual-Chern-Simons Lagrangian. 

In the future we will also consider some additional relevant terms which have not been taken into consideration in 
the present minimal model. As previously shown-2., there generally exists a residual attractive interaction between 
holons and spinons within the t — J model, which should be included when one considers the nodal (d-wave) fermionic 
quasiparticle excitations as "collective" modes in the SC phase. In principle, besides spinon confinement, there also 
exists a holon-spinon confinement in the SC phase of the phase string model, since single spinon or holon excitation 
is not allowed^. How the fermionic nodal quasiparticles can be naturally described in the mutual Chern-Simons 
framework will be a central issue to address in a next study, where the attractive interaction between holons and 
spinons beyond the phase string model should be properly incorporated in order to get a correct excitation spectrum^. 
We do not expect a qualitative modification on the present results of the minimal model by including such a term, 
since the quasiparticles as bound states of holons and spinons are independent, to leading order approximation, of 
those spinon excitations which are confined to form integer neutral spin excitaions discussed in the present work. 
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APPENDIX A: SINGLE- VALUENESS OF $ h (r,r) 

In the transformation l|26() . <5>' l (r, t) is required to be single- valued with mod 2-k in order to ensure the single- 
valueness of the spin operators. Then according to l|27|) . it imposes a constraint on A^, i.e., 

A$ h \c = 2<j> Afap = 2mr, n e Z (Al) 

in which dx M is the tangential differential vector of an arbitrary loop C in (2+1) dimensions. 

Since the gauge field A'^ is an independent dynamic variable, the constraint 1A1|) would be generally violated. 
However, we shall prove below that all the A 1 ^ configurations that violate IjAlf) have vanishing contribution to the 
partition function, which is consistent to the topological constraint © in the Hamiltonian formalism. 

First of all, for an arbitrary loop C in I jAlj l. we may introduce a vortex ring phase configuration e l9 ^ x ^ with an 
arbitrary winding number M (M € Z), which satisfies § D d6 = 2itM for any circuit D that winds around C once, as 
shown in Fig. [21 This singularity in e 10 ^'^ can be clearly expressed by 

e^ x d u d x 0(x) = 2irM j> dy^{C)5^ ( Xfi - y^C)) , M € Z (A2) 

in which y^{C) represents the coordinates on the loop C. 

Then, we can make a singular gauge transformation in terms of such a phase 0(x ll ) as 



h = e iB h, A' = Al + 8^6. 



(A3) 
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FIG. 2: A vortex ring phase 9 is denned such that § d9 = 2nM for any circuit D winding around the loop C once. 



Lagrangian Ch and C s remain invariant, but the mutual-Chern-Simons term in 12Ufl changes as 



Ccs = C cs + -d^ x d v A\ 

7T 



(A4) 



such that the total action is transformed as 



S cS = S eS + d 3 x,-d^ vX d v A\ 



s cS + / d 3 x^-Ay x d v d x e 



7T 

S,. Ll + TIM j A*d.r„. 



(A5) 



Therefore, the partition function can be written as 



Z = J £>[..] exp (-See ~ MM j> A^dx t 

= Const ■ [ D[..} exp ( -S eS - iM I 2A^dx^ 



(A6) 



where D[..] stands for the functional integrations over all the fields, h, h* , z a , z*, A^, and A* The summation over 
M directly lead to the constraint (|A1|) for an arbitrary loop C . 
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